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We investigate the formation of singularities in the incompressible Navier-Stokes equations in 
d > 2 dimensions with a fractional Laplacian |V|". We derive analytically a sufficient but not 
necessary condition for solutions to remain always smooth and show that finite time singularities 
cannot form for a > Oc = 1 + d/2. Moreover, initial singularities become unstable for a > Qc- 
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Scale invariance symmetryjl|, 0, Q holds approximately for the nonlinear Navier-Stokes equations for incompressible 
flow of a Newtonian fluid [1, Q7 The nonlinearity and the scale invariance together can create conditions for the energy 
to cascade down to increasingly finer spatial and temporal scales, e.g., turbulence [5]. In two dimensions, singularities 
cannot form However, more than a century since the discovery of these nonlinear parabolic partial differential 
equations, the question remains unanswered whether or not singularities can form in three dimensions, due to the 
crucial role played by scale invariance. Such fundamental problems remain the subject of ongoing inv estigati ons, due 
to their importance to a number of fields of physics and mathematics 



Here we address the general question of under which conditions dissipation can overcome inertial effects to prevent 
the formation of singularities in finite time. We answer this question by generalizing the problem via a fractional 
Laplacian operator, and then analytically deriving hard inequalities based on the fact that no singularity can form 
provided all partial space derivatives of all orders of the velocity field remain finite for all time. In terms of the 
(i-dimensional Fourier transform v(k,<) of the velocity field v(x,<), we can write 
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To simplify the notation, let {A,B) represent the absolute integral 
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The integrated quantities represent modified Fourier analogs of the statistical moments typically encountered in the 
study of anomalous diffusion and Levy flights [3, [l^, ^21, ^5, 26, 27]. Indeed, our inspiration comes from 

the theory of random walks. In what follows, we use this approach to study the formation of singularities, by deriving 
bounds on the rate of growth of these "absolute moments" ( fc", v ). From Eq. ([T]), we can conclude that no singularity 
can form provided all such absolute moments in Fourier space remain finite. 

The dissipation term in the Navier-Stokes equations contains a Laplacian for incompressible Newtonian fiuids. 
Generalizing the Laplacian operator via Riesz fractional derivatives (2l| . we obtain equations for anomalous 
dissipation JJ: 
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where p denotes pressure, (la the viscosity, / the external forces, p the density, a > 1 the order of the fractional 
Laplacian and i = \ . . .d. Fractional dissipation can arise due to anomalous (i.e., non-Fickian) transport of momentum, 
e.g., in non-Newtonian fluids. Substitution into Eq. ^ of the scale transformation v Av, x A^^x, p X^p, 
and t \T^t, by a factor A, leaves the Navier-Stokes equations unchanged (except for fi). These scale free properties 
leave open the possibility that solutions may conceivably possess structure at arbitrarily small scales, such that for 
a = 2, d > 3 we still do not know whether or not singularities can form in finite time given smooth initial conditions 
and external forces [s], \^ [tI, . Our strategy consists in looking for a spontaneously broken scale invariance symmetry. 
Fourier transforming these equations to eliminate the space derivatives, we obtain 
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with the nonlinearity becoming a nonlocal convolution term in Fourier space. Here, i — \/— T whereas i denotes a com- 
ponent index. We now note that for any function g{t), if dg/dt — A{t) — B(t)g{t) and B > then ^l*?] < |^| — B\g\. 
Identifying g with pvi and B with /Iq/c", we get from Eq. ([5]), 
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and integrating out k over the entire Fourier domain, we get, 
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Using the shorthand notation 
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No approximation has gone into the exact derivation of this hard inequahty. Given smooth initial conditions, the 
foUowing picture thus emerges of the cascade process: excitations of the Fourier modes can pseudorandomly "diffuse" 
anomalously from near the origin k = outwards, towards large k. Wc now ask whether for bounded energy E the 
moments {k",Vi) can diverge in finite time. They cannot diverge at infinite time because nonzero dissipation implies 
that energy decays to zero. 

The terms C{t) do not matter in the context of singularities. We briefly review the known reasons for their 
irrelevance. Due to lack of a separate pressure evolution equation and the special role of pressure in incompressible 
flow, we know that the absolute moments of pressure cannot diverge unless the moments of velocity diverge. Indeed, 
Leray projection methods can eliminate pressure [1,01) because in incompressible flow pressure merely serves to ensure 
Eq. ([3]). It does not contribute to the cascade process, e.g., in d = 3 pressure does not appear in the vorticity equation 
(see below). Similarly, we can ignore smooth external force terms because they do not possess scale invariance. By 
their very definition, smooth external forces have a minimum scale, i.e., they act only within a well defined band of 
spatial frequencies, so ( fc", ) < oo for all n > 0. From now on, we altogether ignore the terms C{t). 

Indeed, the interesting physics only concerns which process wins the competition between the dissipation effect 
that damps high frequency (large k) modes and the nonlocal interaction that transfers (or cascades) the energy to 
the higher frequency modes. In Inequality ([S]), the n-th moment {k"',Vi) can only grow provided the dissipation term 
containing a remains smaller than the inertial terms containing products of moments. We thus arrive at a sufficient 
but not necessary condition for all moments to remain finite: 

for all n,i,j and < £ < n. This condition guarantees that structures of arbitrarily small spatial scale will dissipate. 
However, if the condition does not hold, then Inequality ([5]) may in fact allow a finite time singularity of the form 
(fc"',v) ^ \tc — ^1"'^". The original question of whether or not singularities can form now reduces to the more 
manageable problem of whether or not Eq. ([9]) holds. 

We next look for evidence of spontaneous symmetry breaking. As the excitations in the Fourier modes cascade to 
larger k, the moments grow due to two related but distinct reasons: (i) an "entropic" effect, due to the pseudorandom 
or chaotic spreading or nonlocal diffusion of the excitations in Fourier space and (ii) an energy renormalization effect, 
due to conservation of energy, which leads to approximate normalization of v(k): (v,v) < E/p. The equality 
holds only at t = due to energy dissipation. As an illustrative example of effect (ii), if the energy in a single mode 
with initial amplitude vq becomes evenly split into N"^ modes, then the sum of the amplitudes of the N"^ modes 
grows to N'^/'^vq. Moreover, this energetic growth affects not only (k'^jV) but all moments (fc",v). The interplay 
between these entropic and energetic effects lies at the root of the the fundamental physical mechanism of singularity 
formation. As the moments grow sufficiently large, we know that (fc™,v) ( fc", v ) if m > ti. The highest moment 
appears in the dissipation term, since a > 1. Nevertheless, overgrowth of the numerator of Eq. ([9]) relative to the 
higher moment of order n + a in the denominator may cause the limit not to vanish. We may then lose all control 
and a finite time singularity becomes a possibility. 

The worst case scenario of fastest possible overgrowth for the lower moments in the numerator of Eq. ^ corresponds 
to isotropic cascade with a moment (fc™,v) that diverges in finite time. Moments can diverge only due to fat tails 
(i.e., asymptotic power law decays) in |v(k)|. Hence, we arrive at a multifractal scaling relation for moments in the 
limit of large ( fc™, Vi ) for the worst case scenario: 

n+d/2 

- (fc™,5j)™+'i/2, n>m, (10) 

for all i,j and to > 0. This scaling relation follows directly from conservation of energy and the terms d/2 in the 
exponent come from the energy renormalization effect, as discussed earlier. Substituting Eq. IjlOp into we obtain 
the value of the upper critical ac that guarantees finite moments: 

ac = l + d/2 . (11) 

This critical value constitutes our second important result: singularities cannot form if a > ac. For the marginal case 
a — Oc we can see from Inequality ([8]) that dissipation will hold off the singularity by "buying time" for the energy to 
decay. In fact, the total energy decays always with time, as mentioned earlier. However, we ignore this effect except 
for the marginal case because we can easily show its relative impotence to prevent singularity formation in finite time 
for the non-marginal cases. We thus obtain ac = 5/2 for the important special case d = 3, in agreement with the 
result obtained in refs. 0, [T3 | using different approaches. Eq. pip also correctly gives the known value ac = 2 in 
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d = 2. Moreover, Eqs. ([9]) and (fTO|) together guarantee that pre-existing singularities become unstable under small 
perturbations for a > Uc- 

We next proceed to our third and final result. Critical points often (but not always) separate phases with differ- 
ent symmetries. Intuitively, we expect that a spontaneously broken scale invariance symmetry should lead to the 
emergence of a special characteristic scale that breaks the symmetry. We next show that this indeed happens. If 
a then Inequality ^ shows the existence of a characteristic length scale ijc — 2n/kc for which the inertial 

and dissipative effects can cancel each other. Since only an inequality relation holds, rather than equality relation, 
therefore we cannot make an estimate of this characteristic scale. For a < Uc {a > ac)^ we can make the adimensional 
equivalent of /i^ very large (small) to obtain the dependence of rye on known quantities: 



Note that the units of /i^ depend not only on d but also on a. For a > ac, the length rjc represents the smallest scale 
to which energy can cascade. For a < a^., the value rj^ represents the scale below which inertial effects may dominate, 
conceivably allowing singularities. For d = 3 and a = 2, in the limit of low Reynolds number i?e, we get rjc ~ Re^. No 
matter how viscous or how small we take Re, once a structure forms at length scales smaller than jy^, we cannot rule 
out further cascade to ever finer scales. This somewhat counter-intuitive result suggests that although the parabolic 
Navier-Stokes equations represent a singular perturbation of the hyperbolic Euler equations, yet for certain classes 
of solutions the behavior of ultra-fine structures may possibly remain relatively unperturbed by dissipation at length 
scales 1] ^ r/c for a < Uc- 

We briefly discuss each of these three results. Note that our first result, Eq. ([9]), represents a sufficient but not 
necessary condition for solutions to remain free of singularities. We emphasize the impossibility of obtaining the 
analogous necessary condition from Inequality ([H]). The second result, concerning the upper critical dimension, shows 
precisely the crucial role that dimensionality has on singularity formation. This advance has deep physical significance 
because of the importance of ac in relation to the Hausdorff dimension of the singular set at the time of first blow 
up [3- The third result has a bearing on spontaneous symmetry breaking, which we address next. 

Although the Navier-Stokes equations possess scale invariance symmetry, yet our results, ([Tl]) and (fT2|) . indicate 
spontaneous symmetry breaking. A trivial characteristic length scale ry = always exists for any a, corresponding to 
singular solutions, but our results indicate the existence of a second, nontrivial, length scale r/c. Whereas Inequality 
([8]) does not rule out stable singularities for a < on the other hand for a > ac dissipation always overcomes 
inertial effects at scales smaller than rjc, rendering singularities unstable. The characteristic scale r]c spontaneously 
breaks scale invariance symmetry. In contrast, for a — a^ inertia and dissipation have exactly the same scaling, so 
that the system becomes truly scale invariant, i.e., the original scale invariance symmetry remains unbroken. 

In this context, we note that T. Tao j28i] has pointed out the crucial role of energy and how it behaves under scale 
transformations. If the velocity in a region could evolve to become larger by a factor A, as discussed earlier, then a 
(possibly intermittent) cascade of ever smaller rescaled solutions could repeat the process iteratively until the velocity 
blew up in finite time. For a = 2 and d — 2, the energy remains invariant under scale transformations, hence energy 
considerations by themselves prohibit the formation of singularities. On the other hand, in d = 3 the energy goes 
to zero as A ^ oo, so energy arguments fail. In the general case of Eq. the energy transforms according to 

iJ — > A which we can express in terms of ac &s E X °-i E. Thus, the physical interpretation of ac 

becomes clear: for a > ac singularities become energetically forbidden. 

Finally, we comment on the case a = 2, in terms of the relevant dynamical processes. In d = 3, the vorticity 
cj = Vx V satisfies the vorticity equation 



for inviscid flow [5|. The last term generates vorticity stretching, which can increase vorticity and excite higher 
frequency modes. In d = 2 the vorticity stretching term vanishes, but in d = 3 it allows energy and vorticity to 
cascade and this behavior carries over to viscous flow. The well known criterion of Beale-Kato-Majda fi^ states that 
singularities in the Euler equation in d = 3 can only form if the time integral of the maximum vorticity diverges. 
Whether or not singularities can form in d = 3 for Euler flow, however, still remains an open question. In c? > 3 we 
expect even more violent energy cascades. So Eqs. (ITT|) and suggest that in d = 3 perhaps the more important 
question concerns singularity formation for the Euler equations, since a < ac- Eq. (jl2[) conceivably makes room for 
singularities in Euler flow, if they in fact occur, to persist in the dissipative case. Numerical results have suggested 
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(but not shown) that the Euler equations in d = 3 may in fact allow singularities [l5|. Singularities possibly persist 
even in the Navier-Stokes case [10[ under special circumstances. 

In summary, we have investigated the behavior of solutions of Navier-Stokes equations with fractional Laplacians 
|V|" in d dimensions and shown analytically that singularities cannot form for a > ac- Our approach relies on 
studying the growth of absolute moments (fc", v) in Fourier space. We have used neither mean field approaches nor 
introduced stochasticity. Our results represent hard limits, since they correspond to the worst imaginable scenarios. 
We hope that the approach developed here may inspire future application towards a deeper understanding of other 
nonlinear partial differential equations. 
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